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1 Introduction 



in 

. We study von Neumann entropy and Renyi entropy of spin chains by means 

' of the Fisher-Hartwig formula. The concept of entanglement was introduced 

I Schrodinger in 1935 in the course of developing the famous 'cat paradox' , 

see [68]. Recently it became important as a resource for quantum control, 
which is central for quantum device building, including quantum computers 
. (it is a primary resource for information processing) . Entropy of a subsystem 

^ ' as a measure of entanglement was introduced in [13]. We study spin chains 

with unique ground state. Von Neumann entropy (and Renyi entropy) of the 
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whole ground state is zero, but it is positive for a subsystem [block of spins] . 
In order to define entanglement entropy one has to introduce reduced density 
matrix. The reduced density matrix was first introduced by P. A. M. Dirac 
in 1930, see [27]. 

We calculate the entropy of a block of L continuous spins in the ground 
state of a Hamiltonian. We can think that the ground state is a bipartite 
system \GS) = \AkB), where we call the block by subsystem A and the rest 
of the ground state by subsystem B. The density matrix of the ground state 
is pAB = \GS){GS\, and the density matrix of the block of L neighboring 
spins [subsystem A] is pA = Trg (pab), where we trace out all degrees of 
freedom outside the block. The von Neumann entropy of the block is 

S{pa) = -Tr A {pAln pa), (1) 

which measures how much the block is entangled with the rest of the ground 
state. On the other hand, the Renyi entropy S{pA, oc) is defined as 

5(pA,a) = Y^lnTrA(pS), and a > 0, (2) 

here a is a parameter . Renyi entropy [67] is important in information theory. 
The Renyi entropy turns into von Neumann entropy at a ^ 1. Knowledge 
of the Renyi entropy at arbitrary a permits evaluation of spectrum of the 
density matrix . Our main example is XY spin chain. 

Toeplitz matrix Tl [^] is said to be expressed in terms of generating func- 
tion <?(6') (which is called symbol in mathematical literature): 

rL[f] = (<?i_,), i,j = i,---,L-i (3) 

where ^ 

^k = ^ jj me-'^^'de (4) 

is the fc-th Fourier coefficient of generating function ^(0). The generating 
function ^{9) can be type oi N x N matrix and Tl[5?] is a, N L x N L matrix 
for such case. One of the central objects in the study of Toeplitz matrix Tl [<?] 
is its determinant, which we will denote as £'l[?^], 

Di^[$] := detn[$]. (5) 

Starting with Onsager's celebrated solution of the two-dimensional Ising 
model in the 1940's, Toeplitz determinants play an increasingly central role 
in modern mathematical physics. We refer the reader to the book [60], and 
to survey [59] as for comprehensive sources of the classical results and the 
history concerning the use of Toeplitz determinants in statistical mechanics. 

Another important areas of applications of Toeplitz determinants are ran- 
dom matrices and combinatorics. We refer the readers to the works [71,4,35] 
for the basic results and for the historic reviews. 

Given a generating function ^{0), a principal question is the evaluation 
of the large L behavior of the Toeplitz determinant Di^[<P]. The pioneering 
works on the asymptotic analysis of Toelpitz determinants were done by 
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Szcgo (regular symbol) and by Fisher and Hartwig (singular symbol). These 
results have been used in the study of spin correlation in two-dimensional 
Ising model in the classical works of Wu and McCoy, see for example [60] and 
since then by many other researchers and for a various generating functions. 

The main focus of the majority of works in the area has been, so far, 
the study of spin correlations. The key objects of the analysis have been 
the relevant correlation functions of the local operators. In this paper, we 
discuss yet another, more recent application of the asymptotic analysis of 
Toeplitz determinants in the theory of quantum spin models. Instead of the 
local operators, these applications are concerned with the important nonlo- 
cal objects appearing in spin chains in connection to their suggested use in 
quantum informatics [52]. Indeed, we shall survey some of the recent results 
concerning the quantum, entanglem,ent. We will consider the two applications 
- the entanglement in the XX model and in the XY model. The first one is 
related to a singular scalar generating function, while the second one deals 
with a regular but (2x2) matrix generating function. 

We begin with the brief review of the history and some of the most recent 
results concerning the asymptotic analysis of Toeplitz determinants. 

The plan of the paper is: 

In the second section we discuss the asymptotical expression of the de- 
terminant of a large Toeplitz matrix. The section is divided into subsections. 
Subsection 2.2 is devoted to block Toeplitz determinants. 

Third section is devoted to XY spin chain. In subsection 3.1 we remind 
derivation of determinant representation of entropy of a block of spins in the 
ground state. Isotropic case, i.e. the XX model, is considered in 3.2. For 
anisotropic case we have to use block Toeplitz matrices. 

In section 4 we derive asymptotic expression of entropy of large block of 
spins in isotropic case: the leading logarithmic term and sub-leading correc- 
tions. 

In section 5 we derive asymptotic expression of von Neuman entropy of 
large block of spins in anisotropic case. In the case of XY spin chain the 
entropy has a limit. We calculate the limit. 

In section 6 we calculate limiting expression for Renyi entropy of large 
block of spins in XY spin chain. 

In section 7 we derive the spectrum of the limiting density matrix from 
Renyi entropy. We prove that the spectrum is exact geometric sequence, 
see Eq. 117 and Eq. 124. We also calculate the degeneracy of individual 
eigenvalues, see f Eq. 126. 

The content of sections 4 - 7 is based on the works [45,39,40,33,34]. 

In section 8 we formulate open problems. 



2 Szego and Fisher-HEirtwig Asymptotics 

Throughout the paper we will follow the usual, in the theory and applications 

of Toeplitz determinants, convention to denote the argument of the functions 
on the unite circle either as 6* or as z, z = e*^, i.e. we will always assume the 
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notational identity, 

f{z) = fi0), z = e'0, ee[0,27r). 

We first consider the case of scalar generating function, i.e. N = 1. We shall 
also use for this case the low case symbol instead of 



2.1 Szego and Fisher-Hartwig asymptotics in the case of scalar symbols 

In this subsection we review the basic mathematical facts concerning the 
asymptotics of Toeplitz determinants £>l [0] with scalar generation functions 

Hz)- 

The large L asymptotic behavior of D]^ [0] depends significantly on the an- 
alytical properties of the generating function 4'{6). In the case of the smooth 
enough functions (f){0), the behavior is exponential and its leading and the 
prc-exponcntial terms arc given by the following classical result of Szego, 
known as the strong Szego limit theorem. 

Theorem 1. Suppose that the generation function (f){9) satisfies the condi- 
tions, 

1. (I){e) ^0 , for all e e [0,27r). 

2. index (j>{9) = arg(/)(27r) - arg(/)(0) = 

^- Sfel-oo l^fcl + Sfcl-oo l^ll^feP < where Vk are the Fourier coefficients 
of the function, 

V{9):=\nm, (6) 

that is, 

V{z) = ^fc^"' ^"=2^ I ^i'^y'^^'de. (7) 

k= — (yo ^ 

Then, 



Dl [4>] ^ EsM^^P[LVo), L ^ 00 , (8) 

where the pre- exponential factor, i?sz[<^]; is given by the equation^, 

00 

Es44']=exp(j2kVkV-k)- (9) 
fe=i 

Conditions (1) and (2) on the symbol 4'{9) ensure that the function V{z) 
is a well defined function on the unit circle. Condition (3) is a smoothness 

^ It is this equation which is responsible for the term "strong Szego theorem". 

Szcgo's first result, i.e. Szego limit theorem produced the asymptotics of the deter- 
minant -Dl up to an undetermined multiplicative constant. 
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where 



condition^. It is certainly satisfied by the difiFerentiable functions and is not 
satisfied by the functions having root and jump singularities. In the context 
of Toeplitz matrices, this type of singularities is usually called the Fisher- 
Hartwig singularities. The general form of the symbol ^{z) which has m, 
m = 0, 1, 2, . . . fixed Fisher-Hartwig singularities is given by the equation ^, 

m 

<P{z) = e^(-)^Sr=oft- JJ \z-zj\^"'g.„0,{z)z-^\ z = 9 € [0, 27r), 

3=0 

(10) 

Zj = e*^^ j = 0, . . . , m, = 6*0 < 6»i < • • • < 6»„ < 27r; (11) 

g.M^) ^ g^M = [ < arg. < 2^ ' ^^2) 

^aj > -1/2, /3j e C, i = 0, . . . , m, (13) 

and V{z) is a sufficiently smooth function on the unit circle so that the first 
factor of the right hand side of equation (10) represents the " Szego part" 
of the symbol. The condition on Q.j insures integrability. As it has already 
been mentioned before, a single Fisher-Hartwig singularity at Zj consists of 
a root-type singularity 



\Z-Zi?°'^ 



2sin- 



. 2«i 



(14) 



and a jump gp. (z). A point Zj, j = 1, . . . , m is included in (11) if and only if 
either aj ^ or f3j ^ (or both); in contrast, we always fix zq = 1 even if 
ctQ = Po = (note that gp^iz) = e~^^^°). Observe that for each j = 1, . . . , m, 
z^igf).{z) is continuous at z = 1, and so for each j each "beta" singularity 
produces a jump only at the point Zj. 

In 1968, M. Fisher and R. Hartwig [31] suggested a formula for the lead- 
ing term of the asymptotic behavior for the Toeplitz determinant generated 
by the symbol (10)^. The principal insight of Fisher and Hartwig was the 
observation that the singularities of the symbol yield the appearance of the 



^ In [70], Szego proved this theorem under a somewhat stronger smoothness as- 
sumption on the symbol; namely, he assumed that the symbol is positive, and that 
the symbol and its derivative are Lipshitz functions. It took a substantial period 
of time and the efforts of several very skillful analysts to reduce the smoothness 
conditions to the conditions (1) - (2) above. It also worth noticing that these con- 
ditions are already precise, i.e., if they do not satisfy, the asymptotics (8) might 
not take place. 

In writing the Fisher-Hartwig symbol in form (10) we follow the recent paper 
[22]. Equation (10) is slightly different from the one accepted in most of the liter- 
ature devoted to the Fisher-Hartwig generating functions. The "translation" back 
to the standard form is easy. The main deviation from the standard form is that 
in (10) the product «Sj=o/5j ig factored out which allow to better appreciate the 
non-trivialty of the shifting some of the parameters pj by integers. 

* Some important partial results concerning the asymptotics of the Toeplitz de- 
terminants with singular symbols were also obtained by A. Lenard [53] and used 
by Fisher and Hartwig as a strong evidence in favor of their formula. 
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power-like factors in the asymptotics. Indeed, in the case of all Pj = 0, the 
Fisher-Hartig formula reads as follows. 

Dl M ~ [0]L^r=o "i exp (lFo) , L ^ oo. (15) 

The prc-cxponcntial constant factor, i?pjj[(^], is more elaborated than its 
Szego counterpart £^sz[0] from the Szego equation (8). The description of 
E^pjj [(j>] involves a rather "exotic" special function - the Barnes' G - function 
G{x) which is defined by the equations (see e.g. [75]), 

oo 

G{l + x) = (27r)"/2e-(^+^)^/2-T^^'/2]J{(l + a;/n)"e-^+^'/(2")}, (16) 

n=l 

where 7^ is Euler constant and its numerical value is 0.5772156649 • • •. The G 

- function can be thought of as a "discrete antiderivativc" of the F - function. 
The exact expression for -Efh[0] is given by the equation (cf. Eq. 9), 

fe=l / j=0 

Q<j<k<m 3=0 ^ 

The double product over j < fc is set to 1 if m = 0, so that in the absence of 
singularities, we are back to the strong Szego limit theorem. 

Fisher-Hartwig formula (15) was proven in 1973 by H. Widom [77]. 

The presence of jumps, under the assumption |3?/?j- — 3?/3fc| < 1, does not 
change much the structure of the large L behavior of the Toeplitz determinant 
Di,[0]. Indeed, it is still the combination of the exponential and the power 
terms with the exponential term being determined, as before, by only the 
Szego part of the symbol while the power factor is determined by both the a 
and the (3 parameters of the Fisher-Hartwig part of the symbol. The Fisher- 
Hartwig formula for the general case of symbol (10) reads (cf. Eq. 15), 

£>l[^] ~ Efii[<P]L^^=°^"^'~'^^^ exp(Lyo) , L ^ 00. (18) 

The pre-exponential constant factor, E^fh [</>], is now even more complex than 

in the case of all f3j = 0. In addition to the Barnes' G - fimction, it now 
involves the canonical Wiener-Hopf factorization of the Szego part, e^*^^-*, of 
the symbol (p{z), 

e^(") = b+iz)e^°b-{z), b+{z) = e2^~=i^'="', 6_(z) = e^^"=-~ ^'=^'. 

(19) 

Note that 6+(z) and 6-(z) are analytic inside and outside of the unit circle 
\z\ = 1, respectively, and they satisfy the normalization conditions 5+(0) = 
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6_(oo) = 1. The exact expression for EpHifp] is given by the equation (cf. Eq. 
9 and Eq. 17), 




k=l / j=0 



X 



0<i<fc<m ^ ^ 

~" G(l+a,-+/3,)G(l+a,-/3,-) 



><n ^ c(7;2;o ^'"^ 

The proof of the general Fisher-Hartwig formula (18) is due to E. Basor 

[7] for 3?/3j = 0, E. Basor [8] for a, = 0, |3?/?j| < 1/2, A. Bottcher and B. 
Silbcrmann [17] for \'!kaj\ < 1/2, |3?/?j| < 1/2, T. Ehrhardt [28] for - 
3?/3/c| < 1. The precise statement concerning the large L behavior of the 
Toeplitz determinant -Dl[^] with the Fisher-Hartwig generating function (10) 
is given by the following theorem. 

Theorem 2. (T. Ehrhardt [28]) Let <j){z) be defined in (10), V{z) be C°° on 
the unit circle, ^aj > —1/2, — 3?/3fc| < 1, and aj ± (3j ^ —1, —2, . . . for 
j,k = 0,1, ... ,m. Then, as L ^ oo, the asymptotic behavior of the Toeplitz 
determinant -Dl[(/'] is given by the formulae (18) - (20). 

A. Bottcher and B. Silbcrmann [17] in 1985 and E. Basor and C. Tracy 
[10] in 1991 constructed examples with 3fi/3j not lying in a single interval of 
length less than 1 and such that the large L asymptotics is very different from 
the one given by (18). These examples have showed that for the asymptotics 
(18) to take place, the condition 

\mo-mk\<l, Vi,fc = 0,l,...,m, (21) 

is precise. In the case of arbitrary complex fij, E. Basor and C. Tracy con- 
jectured in [10] a very elegant structure of the large L asymptotics of the 
determinant They based their arguments on the formal analysis of 

the behavior of the both sides of estimate (18) with respect to the shifts 
of the /? - parameters by integers. A detail description of the Basor- Tracy 
conjecture can be found in the original paper [10] as well as in the recent 
work [22] were this conjecture was actually proven with the help of the new 
technique - the Riemann-Hilbert method. 

We refer the reader to monograph [18] and survey [28] for more on the 
mathematics of Toeplitz determinants with the Fisher-Hartwig symbols. 

For the Riemann-Hilbert approach in the theory of Toepitz determinants, 
we refer the reader to the papers [22] and [21] where the method was intro- 
duced (following the similar approach for the Hankel determinants [32] and 
the theory of integrable Fredholm determinants [41]) and to the works [49, 
50,58,57], where the method was further developed. The crucial role in the 
implementation of the Riemann-Hilbert approach to the Toeplitz determi- 
nants is played by the Deift-Zhou nonlinear steepest descent method of the 
asymptotic analysis of the oscillatory matrix Riemann-Hilbert problems [25] 
and by its orthogonal polynomial version [24] . 
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2.2 Block Toeplitz determinants 

A general asymptotic representation of the determinant of a block Toeplitz 
matrix, which generalizes the classical strong Szego theorem to the block 
matrix case, was obtained by Widom in [78,79] (see also more recent work 
[16] and references therein). 

Theorem 4. (H. Widom [79]) Let <?(z) be a N x N matrix function defined 

on the unit circle and satisfying the conditions, 

1. det^(6l) , for all 9 e [0, 2tt). 

2. index det^O) = argdet^>(27r) - argdet<l>(0) = 

3. Er=-ooi^fei+Er=-cx>ifcii^fcp<oo> 

where <l>k are the Fourier coefficients of<P{9), and \F\ denote a, m,a,trix norm of 
the matrix F. Then, the asymptotic behavior of the block Topelitz determinant 
generated by the symbol ${z) is given by the formulae, 

Dl W\ ~ -E'w W\ ( ^ ^ 1^ ^{9)de^ , i ^ oo, (22) 

Ewm=det(Too[^]Too[^-^]). (23) 

where Tao[^] is a semi-infinite Toeplitz matrix, 

Too[^] = (#i-j), i,i = l,2,.... (24) 

From the application point of view, there is an important difference between 
this result and the Szego formula (8) for the case of scalar symbols. Indeed, 
the determinant in the right hand side of Eq. 23 is the Predholm determinant 
of an infinite matrix, and already for 2x2 matrix symbols the question of 
effective evaluation of Widom's pre- factor i?w [^] is a highly nontrivial one, 
even for a relatively simple matrix functions <!>. Indeed, up until very recently 
the only general class of matrix functions for which such effective evaluation 
is possible has been the class of functions with at least one-side truncated 
Fourier series. This class was singled out by Widom himself in [78], and this 
Widom's result has been used in the recent paper [9] of E. Basor and T. 
Ehrhardt devoted to the dimer model. 

Another class of matrix generating functions which admits an explicit 
evaluation of Widom's constant are the algebraic symbols. This fact was 
demonstrated in the works [39, 40, 43] for important cases of the block Toeplitz 
determinants appearing in the analysis of the entanglement entropy in quan- 
tum spin chains. For this class of symbols, Widom's pre-factor admits an 
explicit evaluation in terms of Jacobi and Riemann theta functions. To give 
a flavor of these results, we will now present a detail description of the asymp- 
totics of the block Toeplitz determinant related to the XY spin model ob- 
tained in [39,40]. We shall also use these formulae later in Section 4. 

The Toeplitz determinant in question is generated by the 2x2 matrix 
symbol. 
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where Zi ^ 2:2 are complex nonzero numbers not lying on the unit circle. 
Following the needs of the XY model, we shall assume that the both points 
are from the right half plane though the result we present below can be easily 
generalized to the arbitrary position of the points zi and Z2 outside of the 
unit circle. We will also distinguish three possible locations of the points Zi 
and Z2 on complex plane. 

Case la : Both zi and Z2 are real, they lie outside of the unit circle, and we 
assume that zi > Z2 > I- 

Case lb : Both zi and Z2 are complex, zi = Z2, and we assume that Jf^i > 1 

and 3zi > 0. 

Case 2 : Both zi and Z2 are real, they lie at the different sides of the unit 
circle, and we assume that zi > Z2^ > 1. 

The reason why the Cases la and lb are considered as sub-cases of a single 
case is that in the both these cases all the root singularities of the function 
4>{z) defined in (26) arc inside of the miite circle while all its zeros are outside. 
In Case 2, the zeros and the singularities are evenly distributed between the 
inside and the outside of the unit circle. This difference in the position of 
the roots and singularities of (/>(z) has an impact to the derivations of the 
asymptotics and, as we see below, is reflected in the form of the final answer. 
We shall also see that in the context of the XY model, Case 1 and Case 2 
correspond to the small {h < 2) and large {h > 2) magnetic field, respectively. 

The complex parameter A plays role of a spectral parameter for the 
Toeplitz matrix generated by the symbol, 

Mz)^-^iz)\,^o^[^\^-f^). (27) 

Hence the Toeplitz determinant Dl [^] we are dealing with is in fact a Toeplitz 

characteristic determinant, 

Dl [^] = Dl (A) = det {iXl2L - Tl [#o]) • (28) 

Given the branch points zj of the symbol ^{z), we introduce now the 
elliptic curve, 

w^{z) = (z - zi){z - Z2)iz - z^'Kz - zf 1). (29) 

Let us also re-label the branch points of this curve by the letters A^, A^, 
Ac, and Ad, according to the following rule. Case la : \a = Zi^, — z^^^ 
Ac — -22, Ad = z\\ Case lb : A^ = As = z^^ , Ac = z\, \d = 22, Case 
2 : Xa = z^^, Xb = Z2, Ac = Z2^, Xd = zi. Observe that Xa and Xb are 
always inside the unite circle while Ac and Ad are always outside. This new 
relabeling of the branch points allows to introduce the module parameter of 
elliptic curve (29) in the universal way, 

2 / (30) 



w{zy Jx^ w{z) 
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Theorem 5. ([39,40]) Let 

OC 

6i3(s)= J2 e'^*™'+2'^'"", (31) 



n— — oo 



where t is taken from (30), be the third Jacobi theta- function associated with 
the curve (29). Then, the large L asymptotic behavior of the determinant 
Dl{\) is given by the equations, 

n (X^ ^3 (/3(A) + ^) ^3 (/3(A) -^) , . ^ .... 
Dl(X) ~ -32 (aT\ (1-A) , L^oo (32) 

^3 ( — j 

where 

and a = 1 in Case 1 and a = in Case 2. 

Remeirk. The theta-functions involved in the asymptotic formula Eq. (32) 
has zeros at the points 

±Am, Ato = tanh ( m + ^ tttq, m > 0, (34) 



where, 



TO = -IT = ~i " ; ' > 0. 

r^B dz 



The asymptotics (32) is uniform outside of the arbitrary fixed neighborhoods 

of the points A = ±1 and A = ±A„i. 

Observe that in the case under consideration, det 'P(z) = 1 — A^. Therefore, 
the last factor in (32) is exactly the exponential term of the general Widom- 
Szego formula (22) written for symbol (25). The rest of (32) gives then the 
corresponding Widom's constant, i.e. 

E _ ^3 (/3(A) + ^) e, (/3(A) - ^) 

Similar formulae for the case of the more general quantum spin chains 

were obtained in [43]. The relevant generating function has the same matrix 
structure (25) with the scalar function (/)(z) defined by the equation. 



and p{z) is a polynomial of degree 2n. The analog of the formulae (32) - 
(35) in the case n > 1 involves, instead of elliptic, the hyperelliptic integrals 
and, instead of the Jocobi theta-function, the 2n — 1 dimensional Riemann 
theta-function. 

The methods that lead to these results, involves the theory of integrable 
Fredholm operators [41, 38, 21] and the use of the algebrageometric techniques 
of the soliton theory (see e.g. [12]). 
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3 XY Model and Block Entropy 

The Hamiltonian of XY model can be written as 

oo 

H = - 5^ (1 + 7)<T^<T^+i + (1 - 7)<^^i + ha^„ (37) 

n— — oo 

Here cr^^, are Pauli matrices and /i is a magnetic field; Without loss 

generality, the anisotropy parameter 7 can be taken as < 7; Case with 7 = 
is usually called XX model. The model was solved in [54,1,5,6] and it owns 
a unique ground state \GS). Toeplitz determinants were used for evaluation 
of some correlation functions [6, 69]; Integrable Predholm operators were used 
for calculation of other corrclations[42,26,46]. When the system is in the 
ground state, the entropy for this whole system is zero but the entropy of 
a sub-system can be positive. We calculate the entropy of a sub-system (a 
block of L neighboring spins) which can measure the entanglement between 
this sub-system and the rest part [45]. We treat the whole chain as a binary 
system \GS) = where we denote the block of L neighboring spins by 

sub-system A and the rest part by sub-system B. The density matrix of the 
ground state can be denoted by pab = \GS){GS\. The density matrix of sub- 
system A is PA = TrBipAs)- Von Neumann entropy S{pa) of the sub-system 
A can be represented as following: 

S{pa) = -TrA (pa In pa)- (38) 

This entropy also defines the dimension of the Hilbert space of states of the 
block of L spins. 



3.1 Derivation 

Following Ref. [54] , [55] , wc introduce two Majorana operators 

l-l l-l 
C2i-i = (n<K and C2i = (n<K, (39) 

n=l n=l 

on each site of the spin chain. Operators c„ arc hermitian and obey the anti- 
commutation relations {cm,Cn} = "^Smn- In terms of operators c„, Hamilto- 
nian Hxx can be rewritten as 

N 

Hx.x{h) = i ^(C2„C2„+l - C2„-lC2„+2 + hC2n-lC2n)- (40) 

n=l 

Here different boundary effects can be ignored because we are only interested 
in cascis with N oo. This Hamiltonian can be subsequently diagonalized 
by linearly transforming the operators c„. It has been obtained [54,5] (also 
see [73,52]) that 

{GS\cm\GS) = 0, {GS\c^Cn\GS) = 5 (41) 
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Here matrix Bjv can be written in a block form as 

/ TTo i7_i ... 77i_Ar\ 

ill iio : 



B 



N 



and 



277 70 



d6»e-'"'#o(6'), (42) 



\nN-i Ho ) 

where both i7; and ^oC^*) (for iV 



oo) are 2x2 matrix, 



^o(^) 



m\ 

'-\0) ) 



cos^ — i7sin^ — /i/2 

and (p(9) = — / . ^ — . (43) 

^ |cos^-i7sin6'-/i/2| 



Other correlations such as {GS\cm ■ ■ -CnlGS) are obtainable by Wick theo- 
rem. The Hilbert space of sub-system A can be spanned by rii'^ii'^i'}^' 
where cr^ is Pauli matrix, pi takes value or 1, and vector |0)f denotes 
the ferromagnetic state with all spins up. It's possible to construct a set of 
fermionic operators fej and by defining 



2L 



m 



1, • • • , 2L; bi = {d2i + id2i+i)/2, ? = 1, • • • , L (44) 



with Vmn = (V)„i„. Here matrix V is an orthogonal matrix. It's easy to verify 
that dm is hermitian operator and 

b+ = {d2i - id2i+i)/2, {bi, bj} = 0, {6+ 6+} = 0, {bt,bj} = 5ij. (45) 

In terms of fermionic operators bi and b^ , the Hilbert space can also be 
spanned by ni=i{^i^}^' lOl^jac Here pi takes value or 1, 2L fermionic oper- 
ators bi, bf and vacuum state |0) „ac can be constructed by requiring 

foj|0).ac = 0, ; = 1,---,L. (46) 

We shall choose a specific orthogonal matrix V later. 

Let {f/'/} be a set of orthogonal basis for Hilbert space of any physical 
system. Then the most general form for density matrix of this physical system 
can be written as 

P = ^c(/,J)|V7)(Vj|. (47) 

7, J 

Here c(/, J) are complex coefficients. We can introduce a set of operators 
P{I, J) by P(7, J) a |V'/)(V'j| and J) satisfying 



P{I,J)P{J.K) = Si^k\^ji){'^jiI P{I,J)P{J,K) = Si^K\i^i){-iPi 



(48) 



There is no summation over repeated index in these formula. We shall use 
an explicit summation symbol through the whole paper. Then we can write 
the density matrix as 



p = J2 ^{1, J)P{I, J), c{I, J) = Tr{pP{J, /)). 



(49) 
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Now let us consider quantum spin chain defined in Eq. 37. For the sub-system 
A, the complete set of operators P{I, J) can be generated by OiLi Oi. Here 
wc take operator Oi to be any one of the four operators {b^ ,hi, bfbj, bi bf'\ 
(Remember that bi and b'l are fermionic operators defined in Eq. 44). It's 

easy to find that P(J, I) = (n^=i Oi)^ if P{I^ J) = Ua=i Oi- Here f means 
hermitian conjugation. Therefore, the reduced density matrix for sub-system 
A can be represented as 

PA = J2^rAB (pABillOiA l[Oi. (50) 

\ i=l / i=l 

Here the summation is over all possible different terms Hi^i Oi- For the 
whole system to be in pure state \GS), the density matrix pab is represented 
by \GS){GS\. Then we have the expression for pA as following 

PA = J2iGS\{l[Oi)^\GS)l[Oi. (51) 

This is the expression of density matrix with the coefficients related to multi- 
point correlation functions. These correlation functions are well studied in 
the physics literature [15]. Now let us choose matrix V in Eq. 44 so that the 
set of fermionic basis {bf} and {bi} satisfy an equation 

{GS\bibj\GS) = 0, {GS\bfbj\GS) = Sij{GS\bfbi\GS). (52) 

Then the reduced density matrix pA represented as sum of products in Eq. 51 
can be represented as a product of sums 

L 

PA = '[[(^{GS\b+b,\GS)b+b, + {GS\b,b+\GS)b,b+y (53) 

i=l 

Here we used the equations {GS\b,\GS) = = {GS\b+\GS) and Wick theo- 
rem. This fermionic basis was suggested in Ref. [73,52]. 

Now let us find a matrix V in Eq. 44, which will block-diagonalize cor- 
relation functions of Majorana operators c„. From Eqs. 44 and 42, we have 
the following expression for correlation function of dn operators: 

2L 2L 

{GS\dmdn\GS) = J2J2''rni{GS\CiCj\GS)vjn , 

i=l j=l 

{GS\ \GS) = 5jnn + i(BL)mn, 

{GS\dMGS) = (5„„ + i(BL)™„. (54) 

The last equation is the definition of a matrix Bl. Matrix Bl is the sub- 
matrix of Bn defined in Eq. 42 with m, n = 1, 2, . . . , L. We also require Bl 
to be the form [73,52] 

Bl = FBlF^ = eL=ii^n^ f -1 J ) = " ® f -1 J ) • (55) 
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Here matrix is a diagonal matrix with elements Vm (all Vm are real num- 
bers). Therefore, choosing matrix V satisfying Eq. 55 in Eq. 44, we obtain 
2L operators {bi} and {b'^} with following expectation values 

{GS\bm\GS) = , {GS\bmbn\GS) = , (C?5|6+6„|G5) = (56) 

Using the simple expression for reduced density matrix pA in Eq. 53, we 
obtain 

PA = 11 (^btb, + '-^hbA . (57) 

i=i ^ ^ 

This form immediately gives us all the eigenvalues XxiX2---xi, of reduced den- 
sity matrix pA, 

Axix....x. = n Xi = 0,iyi. (58) 

i=l 

Note that in total we have 2^ eigenvalues. Hence, the entropy of pA from 
Eq. 38 becomes 

L 

5(p^) = ^e(l,z.„) (59) 

m=l 

with 

x + u, ,x + v. x — v^ ,x — v. 
e{x,u) = _ln(^-)-^-ln(^-). (60) 



3.2 XX model 

Notice further that for XX model, i.e. 7 = case, matrix Bl can have a 
direct product form 



Bl = Gt 



1 
-1 



with Gl 



/ 00 0- 



\(t>L-l ... 



(61) 



where (bi is defined as 



2lT 



1, -kp <0 <kF, 
-1, kp <0 < {2iT-kF) 



(62) 



and kp = arccos(|/;,|/2). From Eqs. 55 and 61, we conclude that all Vrn are 
just the eigenvalues of real symmetric matrix Gl- 

However, to obtain all eigenvalues Um directly from matrix Gl is a non- 
trivial task. Let us introduce 



I)l(A) = det(GL(A) = A/l - Gl) . 



(63) 
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Here Gl is a Toeplitz matrix and /l is the identity matrix of dimension L. 
Obviously we also have 

L 

DL{X)=l[{X-iym). (64) 

m=l 

Prom the Cauchy residue theorem and analytical property of e{x,i'), then 
S{pa) can be rewritten as 

S{pa) = lim / dA e(l + e, A) In (A) . (65) 

Here the contour F' in Fig 1 encircles all zeros of -Dl(A) and function 
e(l + e, A) is analytic within the contour. Just like Toeplitz matrix Gl is 




Fig. 1 Contours F' (smaller one) and F (larger one). Bold lines (— oo,— 1 — e) 
and (1 + e, oo) are the cuts of integrand e(l + e, A). Zeros of Dl{X) (Eq. 64) are 
located on bold line (—1,1). The arrow is the direction of the route of integral we 
take and r and R are the radius of circles. ^ 



generated by function (/)(0) in Eqs. 61 and 62, Toeplitz matrix Gl(A) is gen- 
erated by function ^{6) defined by 

'f'^^^ = \X + l,kF<e<{2n-kF). 

Notice that 4>{^) is a piecewise constant function of on the unit circle, 
with jumps at ^ = zbfci?- Hence, if one can obtain the determinant of this 
Toeplitz matrix analytically, one will be able to get a closed analytical result 
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for S{pa) which is our new result. Now, the calculation of S{pa) reduces to 
the calculation of the determinant of Toeplitz matrix Gl(A). 



3.3 XY model 

Similarly let us introduce: 

Bi(A) = i\lL - Bi, Dl{\) = det Bi(A). (67) 
Here II is the identity matrix of dimension 2L. By definition, we have 

DL{X) = i-l)^f[iX'-t.l). (68) 

m=l 

Using again the Cauchy residue theorem we obtain that, similar to (65), 

S{pa) = lim -ir / dAe(l + e,A)-^lnI?i(A) . (69) 
e^o+ 47n Jp' dA 

Here the contour F' in Fig 1 encircles all zeros of Dl(X). 

We also realized that Bz,(A) is a block Toeplitz matrix with the generator 

${z), i.e. 

/ no n-i ...ni-L\ 



Bl(A) 



ill ilo 



with 



and (l){z) 



\nL-i iio / 

dzz-^-'^^{z), ${z): 

XI (l-Aiz)(l-A2Z-i) 



iA <p{z] 
-(p-Hz) iA 

1/2 



(70) 



(71) 



(72) 



^Ai (l-Af^-i)(l-A^ z)^ 

We fix the branch by requiring that (f>{oo) > 0. We use * to denote 
complex conjugation and S the unit circle shown in Fig. 2. Ai and A2 are 
defined differently for different values of 7 and h. There are following three 

different cases: 

In Case la {2^/1 — 7^ < h < 2) and Case 2 {h > 2), both Ai and A2 are 

ft- Vfe2_4(l_^2) 1 + ^ 
'^l = , ^ , ^^2 = M- 



2(1 + 7) 1-7 
In Case lb (ft^ < 4(1 — 7^)), both Ai and A2 are complex 

ft - 1^4(1 - 72) - /i2 



Ai 



2(1 + 7) 



A2 = 1/At. 



(73) 



(74) 



Note that in the Case 1 the poles of function (f){z) (Eq. 72) coincide with the 
points Xa and As, while in the Case 2 they coincide with the points A^ and 
Ac. 
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z-plane with z = x + i y 




Case la&2 Case lb 



Fig. 2 Polygonal line E (direction as labeled) separates the complex z plane into 
the two parts: the part fi+ which lies to the left of E, and the part O- which lies 
to the right of S. Curve E is the unit circle in anti-clockwise direction. Cuts Ji, J2 
for functions cj>[z),w(z) are labeled by bold on line E. Defimtion of the end points 
of the cuts A... depends on the case: Case la: Xa = Ai and As = AJ^, Ac = A2 
and \d ~ A^^. Case lb: Xa ~ Ai and Xb ~ AJ^, Ac ~ X^^ and Xd = A2. Case 
2; Xa = Ai and Xb ~ X2, Ac ~ A^ and Xd ~ X^ . ^ 



4 Block entropy of XX model and Fisher-Hartwig Formula 

From Eq. 65, one needs the calculation of Toeplitz determinantZ?i(A) with 
a singular generating function 

A - 1, -kp <0 < kp, , , 

X+l,kF<0<{27T-kF). ^ ' 

It is easy to check that this function admits the canonical Fisher-Hartwig 
factorization given by Eq. 10 with 

TO = 2, a,^0 Vj, /3o = 0, /32 = = /?(A) - ^ In (76) 

Zm A — 1 



A + 1^ -'''^^^ 



and 

e^(-)^e^o = (A + l)(^^j . (77) 

The branch of the logarithm is fixed by the condition, 

< arg (x^) < 

For A ^ [—1, 1], the left inequality is also strict, and hence |3fi(/3i(A))| < ^ and 
|5R(/32(A))| < ^. Therefore, Theorem 2 is applicable (indeed, even its earlier 
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weaker version proven by E. Basor [8] would be suffice) and we see that the 
determinant -Dl(A) of A/l — Gl can be asymptotically represented as 



L'l(A) = (2-2cos(2M) " ^^\g(i + f3{\))G(l- (3{\))Y 

(A + 1)((A + 1)/(A-1))"''"^"| (79) 

Here G is, as before, the Barnes G-function and 

G(l + /3(A))G(1 - /3(A)) = e-(^+^-)P'W ft { " ^) " e'^''^'^"' } • 

(80) 

Let us substitute the asymptotic form Eq. 79 into Eq. 65 and after some 
simplification [45], we have that 

5(p^) = ilnL+iln||l-Qyj+^+ri, L ^ oo (81) 

with 

io I 3i tsinh2(t/2) 2sinh3(V2)J ^ ' 

for XX model. The leading term of asymptotic of the entropy ^ InL in Eq. 81 
was first obtained based on numerical calculation and a simple conformal 
argument in Ref. [73,52] in the context of entanglement. We also want to 
mention that a complete conformal derivation for this entropy was ftnmd 
in Ref. [48]. One can numerically evaluate Ti to very high accuracy to be 
0.4950179 • • •. For zero magnetic field {h = 0) case, the costant term Ti + 
In 2/3 for S{pa) is close to but different from (7r/3) ln2, which can be found 
by taking numerical accuracy to be more than five digits. 



{< 



5 Block entropy of XY model and block Toeplitz determiniant 

For the block entropy of XY model, by virtue of Eq. 69, our objective be- 
comes the asymptotic calculation of the determinant of block Toeplitz matrix 
Dl{X) or, rather, its A -derivative ^lnZ)i^(A). 
Let us denote, 

zi := Xi^, and Z2 := A2. (83) 

It is easy to check than that the generating function introduced in Eq. 71) - 
Eq. 72 coincides with the one introduced in Eq. 25) - Eq. 26 together with the 

case-separations and the Xa - Xd labeling of the branch points. Hence one 
can use Theorem 5 and substitute the asymptotic form Eq. 32 into Eq. 69. 
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Deforming the original contour of integration to the contour F as indicated 
in Fig. 1 we arrive at the following expression for the entropy [39,40]: 

wich can also be written in the form, 

'^''^-2i WW) J^ii^ ^''^ 

This is a limiting expression as L ^ oo. In [40] it is also proven that the 

corrections in Eq. 84 are of order of O ^A^^/Vl^ . 

The entropy has singularities at phase transitions. When r ^ we can 
use Landen transform (see [75]) to get the following estimate of the theta- 
function for small r and pure imaginary s: 



^ taT\ = —s T Triers + O 5— s , as r ^ 0. 

Now the leading term in the expression for the entropy (84) can be replaced 

by 

iTT _ /e"*'^/^ 



S{pa) = -^+0 [—^j ^ ^ 0- (86) 

Let us consider two physical situations corresponding to small r depending 

on the case defined on the page 2: 

1. Critical magnetic fi,eld: 7 7^ and h ^ 2. 



This is included in our Case la and Case 2, when h > 2^ \ — ^'^ . As h ^ 2 
the end points of the cuts \b Ac, so r given by Eq. (30) simplifies and 
we obtain from Eq. (86) that the entropy is: 

S'(p^) = -lln|2- /i| + ^ln47, for /i ^ 2 and 7 ^ (87) 

correction is 0(|2 — /i| In^ |2 — ft|). This limit agrees with predictions of 
conformal approach [48,19]. The first term in the right hand side of (87) 
can be represented as (1/6) In ^, this confirms a conjecture of [19]. The 
correlation length ^ was evaluated in [5] . 
2. An approach to XX model: 7 — > and /i < 2: It is included in Case lb, 
when < /?- < 2-y/l — 7^. Now Ag ^ Ac and A^ ^ A^, we can calculate 
T explicitly. The entropy becomes: 

S'0(p^) = _lln7+lln(4-/i2) + ^ln2, for 7 ^ and /i< 2 (88) 
correction is 0(7 In^ 7). This agrees with [45] (see also Eq. 81). 
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As it has already been indicated, the theta-functions involved in the 
asymptotic formula Eq. (32) has zeros at the points ±A„j which are defined 
in Eq. (34). Theorem 5 shows, in particular, that in the large L limit, the 
points ±Ato are double zeros of the Dl{\). More precisely, we see that in the 
large L limit the eigenvalues and fc'2m+i from (59) merge to A^: 

l'2mj 2^2m+l — * Am, (89) 

which in turn implies (cf. Eq. 59) the following equivalent description of the 
limiting entropy S{pa) [39]: 

The limiting entropy, S{pa), of the subsystem can he identified with the 
infinite convergent series, 

oc oo ^ 

S{pa)= J2 e(l'^™)= E + ^ (90) 

m=—oo m= — oo 

Indeed, equation (90) follows from the substitution of Eq. (32) into Eq. (69) 
and integrating over the original contour F' of Fig. 1 

It is also worth mentioning that relation (89) also indicates the degeneracy 
of the spectrum of the matrix Bj, and an appearance of an extra symmetry 
in the large L limit. 

Remeirk. These numbers A^ satisfy an estimate: 

\Xm+i - Ami < 47rro with tq = -ir. 

This means that (A^+i — Am) — > as r — > for every m. This is useful for 

understanding of large L limit of the XX case corresponding to 7 0, as 
considered in [45] . The estimate explains why in the XX case the singularities 
of the logarithmic derivative of the Toeplitz determinant dlnDL(X)/dX form 
a cut along the interval [—1, 1], while in the XY case it has a discrete set of 
poles at points ±Xm of Eq. (34). 

The higher genus analog of formula Eq. (84) for the class of quantum spin 
chains introduced by J. Keating and F. Mezzadri in [47] has been obtained 
in [43]. 

Remark. It was shown by Pcschcl in [65] (who also suggested an alternative 
heuristic derivation of equation (90) based on the work [19]), the series (90) 
can be summed up to an elementary function of the complete elliptic integrals 
corresponding to the modular parameter t - see Eqs 109 and 110 below. It is 
an open problem whether an analogous representation of the integral Eq. (84) 
exists for higher genus. The key issue here is the extreme complexity of the 
identification of the zero divisor of the theta- functions in the dimension 
grater than 1. 

6 Renyi entropy and the spectrum of reduced density matrix of 
XY model 

The Renyi entropy of Sa{pA) of the block of spins is defined by the expression 

SJpa) = InTrip'X), a^l and a > 0. (91) 

1 — a 
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Here the power a is a parameter. The Renyi entropy is intimately related to 
the spectrum of the reduced density matrix pA- Indeed, let A„, (0 < A„ < 1) 
and a„ denote the eigenvalues and their multiplicities of the operator pA. 
The spectrum is completely determined by its momentum function, i.e. by 
the C-function of pA , 



The obvious equation takes place, 

CpA (a) 



n=0 



_ Jl-a)SRipA,a) 



(92) 



(93) 



The key point is that we can evaluate Sa{pA), and hence CpAi'^)^ explicitly. 

As it is shown in [45], the Renyi entropy So:{pa) of a block of L neighbor- 
ing spins, before the large L limit is taken, can be represented by the finite 

sum, 

L 

/ I -1- 7y. \ /I — 7/7 \ 

(94) 



1 ^ 

SnipA, a) = V In 

1 — a f—f 



where the numbers 



1,...,L 



are the eigenvalues of the same block Toeplitz matrix Eq. 68 as we worked 
with in Section 3.4. In virtue of Eq. (89), the Renyi entropy in the large L 
limit can be identified with the convergent series. 



SR{pA,a) 



1-a 



E 1- 



1 + A„ 



(95) 



with 



Xm = tanh ym + — ^— ) '^to- 



(96) 



The summation of the series can be done following the same approach as in 
[65] in the case of the von Neuman entropy. The result is (for details see [33]) 
the following. 



SR{pA,a) = 



1 - a V 12 6^4 



-aTTTo 



^lnn(l + e+^)', (97) 

(98) 



for the case h> 2, and 

Sr{pa,oi) = -r^— 
1 — a 



^ro ^ 1 , k'\ 



n=0 



oo 



n=l 



\ — a 



In 2, 



(99) 
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for the case h < 2. In these equations, tq = —ir is the module parameter 
defined in Eq. (30,) and k = k(qi), k' = k'{qi) arc the standard cUiptic 
modular functions, see e.g. [75]. The quantities k and k' are simply related 
to the basic physical parameters 7 and h. Indeed, one has that 



^(V2)2+72— /^^ Case la: 4(1 - 7^) < /i^ < 4; 

k = \ y^(l -/^V4-72)/(l -/t74) , Case lb: /i^ < 4(1 - 7^); (100) 
, 7 / ^{hl2y + 72-1 , Case 2 : /i > 2. 

k' = \J\ - fc2. 

By standard techniques of the theory of elliptic functions, equation Eq. (30) 
can be transformed into the following representation for the module tq as a 
function of k. 



To 



I{k') 



k' = Vl - fc2, (101) 



I{k) ' 

I{k) is the complete elliptic integral of the first kind, 

I{k)=C^ ''^ (102) 

The g-products in Eq. 97 and Eq. 99 can be expressed in terms of the 
elliptic lambda function or A - modular function. The A - function plays a 
central role in the theory of modular functions and modular forms, and it is 
defined by the equation (see e.g. [75]), 



ir) 



Mr) = ^Dt^ = fc'(e'""), Qr>0, (103) 



where 9j{s\T), j = 3,4 are Jacobi theta-fucntions; the function ^3(s|t) has 
already been defined in Eq. 31, while the function ^4(s|t) is defined by the 
equation, 

oc 

6»(s|r) = (-l)"e''*™'+2''^'*". (104) 

n= — 00 

The A - function is analytic function of r, Sjt > 0, and it satisfies the follow- 
ing symmetry relations with respect to the actions of the generators of the 
modular group, 

A(-i) =1-A{t). (106) 

In terms of the A - modular function, the formulae for Renyi read as 
follows [33]. 

SR{pA,a) (107) 

= I ln(fc k') - ^ ln(A(mro)(l - A(mro))) + iln2, 

o 1 — a 12 1 — av / 6 
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for /i > 2 and 



Sr{pa, a 
1 a 



(108) 



6 1-a 



In 



In 



A^(Q!zro) 



12 1-a l-A(mTo) 3 



In 2, 



for /i< 2. 

Eqs. 107 and 108 allow to apply to the study of the Renyi entropy the 
apparatus of the theory of modular functions. 

Remark. Using Eq. 107 and Eq. 108 one can evaluate the asymptotics of 
the Renyi entropy as a ^ 1. This would lead to the following formulae for 
the Neumann entropy, 



S{pa) = 



1 



ln( ) + (1 



in Case 1, and 



S{pa) = 



12 



ln( 



P 4I{k)I{k') 

2 ' TT 



In 2, 
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(fc2fc'2; 



+ {k' - k'^) 



aAI{k)I{k') 



(109) 



(110) 



in Case 2. For the Cases la and 2 these formulae were first obtained by 
Peschel in [65] by a direct summation of series (90) 



7 Spectrum of the limiting density matrix 

We will show now, following [34] , how to extract from Eq. (97) and (99) the 
information about the spectrum of the density matrix pA- 

Consider first the case h > 2. Combining equations (97) and (93), we 
arrive at the following representation for the C- function ^^^(a), 

Cp.(a)=e"(^ + ^'""^)n(l+e^^)'- (111) 

n=0 

At the same time, using the classical arguments of the theory of partitions 
(see e.g. [2], Chapter 11, equation (11.1.4)) we have that 

oo oo 

l[{l+q'-+')^J2Po\n)q\ (112) 

n=0 n=l 

where ^^^(0) = 1 and pQ^(n), for n> 1, denote the number of partitions of 
n into distinct positive odd integers, i.e. 

#{(mi,...,mfe) : mj = 2rj + 1, toi > TO2 > ... > ruk, 

n = TOi + 7712 + ... + mfe} . 

Hence (111) becomes, 

, / ^ oo 

C,,(a) = eH^+^i"^J ^a„g2, (113) 

n=0 



24 



where, 

n 

ao = 1, an = Y.Po\l)Poin - 0- (114) 

1=0 

Finally, observing that 

q2 = {e-^rony^ (115) 

we conclude that 

oo 

Cp,(a) = ^a„A^:, A„ = e-'^-''"+^+^^"^. (116) 

Comparing the last equation with equation (92) we arrive at the following 
theorem. 

Theorem 6. ([34]) Let the magnetic field h > 2. Then, the eigenvalues of 
the reduced density matrix pA are given by the equation, 

A„ = e-'^^°"+'S^+5'"^, n = 0,1,2,...., (117) 

and the corresponding multiplicities a„ are determined by the relation (Hi). 

The case ft < 2 is treated in a very similar way. Instead of (111) we have 
now the formula, 

C,.(a)=2eK-^+^'"A.)n(l + e)', (118) 
where Qa as in (98). The analog of the Taylor expansion (112) is the equation, 

oo oc 

n(i+«^")=E^'^^w9'"' (119) 

n— n=l 

where p^^(O) = 1 and Pj^\n), for n > 1, denote the number of partitions of 
n into distinct positive integers, i.e. 

# {(mi, ...,mk) : mi > m2 > ... > mfc > 0, n = mi + m2 + ... + m^} . 

Hence (118) becomes. 



where. 



CM = 2e"V-^-^-i^; (120) 

n=0 



bo = l, 6„ = Ep^)(0p^^(n-0. (121) 
z=o 

Finally, observing that 



= (e-2--on')"^ (122) 
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we conclude that 

oo 

U(a) = 2^6„A«, A„ = e-2'^^<'"-^+^^"^. (123) 

Comparing the last equation again with Eq. 92 we arrive at the analog of 
Theorem 6 for the case h < 2 . 

Theorem 7. ([34]) Let the magnetic field h < 2. Then, the eigenvalues of 
the reduced density matrix pA are given by the equation, 

An = e-2--o"-^+|i"3li, n = 0,l,2,...., (124) 

and the corresponding multiplicities equal 26„ where the integers bn are de- 
termined by the relation (121). 
Let 

OO 

/(ar) :=5^a„a;", (125) 

be the generating function for the coefhcients a„. Then, Eq. 107 and Eq. 
93 in conjunction with the symmetry property (106) allow to analyze the 
asymptotic behavior of the functio f{x) generating function as a: ^ 1. In 
its turn, this fact yields the evaluation of the large n asymptotics of the 

multiplicities a„ (details are in [34]). 

Theorem 8. Let an be the multiplicities of the eigenvalues of the reduced 
density matrix for h> 2. Then their large n behavior is given by the relation, 

On ~ 2-3/23-i/4n-3/Vv^, n ^ OO. (126) 



8 Summary and Open Problems 

We want to emphasise that the method described here also works for eval- 
uation of correlation functions. For example space, time and temperature 
dependent correlation function of quantum spins was evaluated in [42] . The 
book [15] explains how to apply this method for calculation of correlation 
functions in Bose gas with delta interaction. 

On the other hand there are still open problems. For example let us 
consider the XX Z model. The Hamiltonian can be written in terms of Pauli 
matrices (T„: 

oo 

Hxxz = - (127) 

n— — oo 

At Z\ < — 1 the model has a gap and the ground state is anti- ferromagnetic. 
Challenging problem is to calculate the von Neumann entropy and Renyi 
entropy of large block of spins on the infinite lattice. It will be interesting to 
find the dependence of limiting entropy on A. 
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